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ABSTRACT 

The utility of the ribosomal DNA gene complex for species identification of Actiniaria was 
examined. The use of universal ribosomal PCR primers is problematic in this group due to 
the presence of algal symbionts. Universal primers were initially used to amplify a region 
containing partial 18S, complete ITS, 5.8S, ITS2, and partial 28S sequences from six 
sea anemone species. The development of two sea anemone specific primers for this 
region was necessary to avoid amplification of algal symbionts for a number of species. 
Complete sequences of the 18S-28S fragment were obtained from three species, 
?Anemonia sp. (724 bp), Heteractis malu (670 bp) and Stichodactyla haddoni (734 bp); 
partial or non-overlapping sequences were obtained from Entacmaea quadricolor 
(480bp from 18S), Macrodactyla doreens/s (523 bp: 300bp from 18S and 223bp from 
28S) and Oulactis muscosa (556bp: 285bp from 18s and 271bp from 28S). Average 
sequence divergence among sea anemone species was approx. 24% indicating that this 
region may indeed be useful for species identification. However, unexpectedly low 
divergence recorded between two species in different genera, neither of which could be 
verified by histology due to specimen unavailability, indicated that traditional histological 
methods are still needed to confirm identification and certainly until such time that an 
rDNA database of sea anemone tissue has been established. □ ribosomal DNA; sea 
anemone specific primers; universal primers 


In February 2005, the Australian Marine 
Sciences Association, SEQ Branch, hosted the 
T hirteenth International Marine Biological Work¬ 
shop, The Marine Fauna and Flora of Moreton 
Bay, Queensland. Fieldwork was conducted over 
a period of three weeks and occurred in a variety 
of environments including off-shore reefs, small 
islands accessible only at low tide, piers, estua¬ 
rine mouths and mud flats. A taxonomic paper 
documenting the species found is presented by 
Fautin el al. (2008, this volume). Of the more 
than 20 species that are now known from More¬ 


ton Bay, we obtained tissue from the following six 
species and genera of anemones to assess the 
usefulness of DNA in identification, and the poten¬ 
tial for understanding phylogenetic relationships: 
lAncmonia sp., Heteractis main, Stichodactyla 
haddoni, Entacmaea quadricolor, Macrodactyla 
dorcensis and Oulactis muscosa. 

Species identification of sea anemones (Antho- 
zoa: Actiniaria) can be difficult, especially in the 
field. The taxonomic key currently utilised, de¬ 
signed by Oskar Carlgren (1949), is based mainly on 
histological differences and therefore requires 
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Table 1. Moreton Bay sea anemone species used in the genetic analysis. 


Collection 

Location 

Latitude 

Longitude 

Field 

Identification 

Laboratory 

Identification 

Genetic Source 

Museum 

Reg. No. 

Bird Island 

27° 30' S 

153° 23' E 

Eutacmaea 
quadricolor juv. 

?Anemonia sp. 

Whole animal 

-s 

MTQ G58754 

Flat Rock, N. 
Stradbroke I. 

27" 24' S 

153° 33' E 

Entaamca 

quadricolor 

Not available 

Tentacle 

Whole spec, 
not coll. 

Dunwich, flats in 
front of MBRS 

27" 30'S 

153° 24'E 

Stichodactyla 

haddoni 

Not available 

Pedal disc & 
tentacle sample 

Whole spec, 
not coll. 

Shag Rock, N. 
Stradbroke I. 

27° 24.85'S 
153°31.59'E 

Heteractis main 

Heteractis 

main 

Pedal disc & 
tentacle sample 

MTQ G58749 
QM Unreg. 

Frenchmen's 

Beach 

27° 25' S 

153° 32' E 

Oulactis 

muscosa 

Oulactis 

muscosa 

Pedal disc 

MTQ G58756 

Dunwich 

27" 30' S 

153° 24' E 

Macrodactyla 

doreensis 

Macrodactyla 

doreensis 

Pedal disc & 
tentacle 

MTQ G58748 


collection of whole animals, which may not always 
be practical. Furthermore, histological analysis of 
sea anemones is time consuming and requires 
considerable expertise as some closely related 
species are almost impossible for the non specialist 
ntify, often resulting in incorrect taxonomic 
assignment (Stephenson 1928; Fautin 2000; Hausser- 
mann 2004). 

Identification is further complicated by the fact 
that some species are virtually identical in appear¬ 
ance, distinguished by only one or two morpho¬ 
logical features. For example, the two species 
Heteractis main and Heteractis crispa are differ¬ 
entiated in the field on the basis of 1) column 
texture, which is firmer (leathery) in H. crispa 
than in H. main, and 2) tentacle length, which is 
meant to be twice as long in H. crispa than H. main 
(Fautin & Allen 1997). Both these characteristics 
can be misleading since the former is open to 
subjective interpretation if both species arc not 
present side by side in the wild, while the latter 
may not necessarily be useful as the tentacles can 
be contracted at the time of collection/observation. 

Alternatively, delineation of some species may 
be quickly achieved using appropriate molecular 
genetic methods (eg the Barcoding of Life initiative 
Hebert cl al. 2003). However, previous genetic 
studies including sea anemone taxa have either 
only focused on questions pertaining to higher 
order anthozoan relationships (Won et al. 2001, 
Daly et al. 2003) or intraspecific population structure 
(Hunt & Ayre 1989). Numerous mitochondrial DNA 


genes such as COl (Fautin & Smith 1997), 
COI1I and 16S rDNA (Geller & Walton 2001) 
have been used to infer phylogenetic relation¬ 
ships among the Actiniaria. However, mitochon¬ 
drial gene sequence divergences within and 
among anthozoan families, including sea 
anemones, has been found to be significantly 
lower than other marine invertebrate species 
(Shearer et al. 2002). Barcoding studies also 
discovered that mitochondnal DNA evolved too 
slowly in sea anemones and other cnidarians 
for mtDNA differences to be an informative 
indicator of species (Hebert et al. 2003). Inter¬ 
estingly, Shearer et al. (2002) also found that, 
unlike all other metazoan taxa, substitution 
rates in anthozoan nuclear genes are much 
higher than in mitochondrial genes and there¬ 
fore may be of greater utility in terms of species 
identification. Indeed, a number of other studies 
have suggested that the nuclear ribosomal 
(rDNA) gene complex incorporating 18S, ITS1, 
5.8S, 1TS2 and 28S could be ideally suited to 
examining below genus level relationships within 
the Actiniidae (McCommas 1991; Odorico & 
Miller 1997). Most recently Acuna et al. (2007) 
used the ITS region of rDNA in addition to 
morphology to distinguish between tliree species 
within the genus Aulactinia. 

Molecular studies of sea anemones can be 
potentially complicated by the presence of sym¬ 
biotic algae or zooxanthellae in the anemone 
tissue (Shearer et al. 2005) and possibly tissue 
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consistency (Pinto et ai 2000). If species possess 
zooxanthellae they are generally found in the 
gastrodermal tissues (i.e. tentacles and oral discs), 
although in some species they can be distributed 
heterogeneously throughout their hosts, being 
rare in only the pedal disc region or mesenteric 
tissue layers (Fautin & Smith 1997; Hciusser- 
mann 2004). Therefore DNA extractions can con¬ 
tain both the host and algal genomes, which 
may cause confounding results especially for 
sequence data generated using broadly con¬ 
served or 'universal' primers (see Shearer et ai 
2005). A study by Pinto etol (2000) found tissue 
consistency to impinge on the success of extrac¬ 
tion of DNA from sea anemones, due to hard¬ 
ness of tissue from being preserved in ethanol. 
They concluded that a slow and gradual diges¬ 
tion method was optimal for extraction. 

Here we conduct a preliminary study to 
examine the utility of the rDNA gene complex 
in the identification of sea anemone species and 
test whether a known universal primer pair is 
sufficient for such studies or whether anemone 
specific primers will be required. Furthermore 
we use modern DNA extraction kits to see if 
previous problems associated with sea anemone 
DNA extraction can be circumvented. 

MATERIALS AND METHODS 
SPECIMEN AND TISSUE COLLECTION 

Collection techniques included; removing 
anemones from rocks by chisel and hammer, 
scraping animals off rocks by fingernail or taking 
a small tissue sample from the animal in the 
wild for genetic analysis if identification was 
100% positive in the field. Tissues for analysis 
were collected from twelve species (based solely 
on field identifications). Of these, six samples 
representing an initial five species were used in 
the genetic analysis (Table 1). Additional samples 
of Heteractis mcihi were collected from Shag Rock 
subsequent to the workshop. 

In order to examine and minimise possible 
zooxanthellae contamination, small tissue samples 
of less than 5mm in length were excised from 
either the lower column/pedal disc or, where 
possible, separate tissue samples from both the 
tentacles and pedal disc region of each species 
were taken. All samples for genetic analysis 
were stored in 100% ethanol. Where whole 


specimens were collected, tissue samples were 
taken after animals were relaxed in magnesium 
chloride and before being preserved in 10% 
formalin: seawater. All ethanol preserved tissue 
samples were stored at -20°C until genetic anal¬ 
yses were performed. Hetevactis wain specimens 
collected subsequent to the workshop were stored 
in 100% ethanol and kept at room temperature 
(approximately 21 °C) only. 

DNA EXTRACTION, PCR AND SEQUENCING 

To test to the usefulness of modern DNA 
extraction kits with ethanol preserved sea anem¬ 
one tissues, total genomic DNA was extracted 
from both tentacles and pedal disc tissues using 
DNeasy Tissue Kits (QIACEN) as opposed to 
the far more labour intensive protocol of Pinto 
et ai (2000). Partial 18S rDNA, complete ITS1, 
5.8S, ITS2 and partial 28S rDNA sequences were 
initially amplified using the primer pairs RA2 
and 1TS2.2 described by Wbrheide (1998) RA2 
is located in the flanking 3' end of the small 
subunit ribosomal gene (18S) and ITS2.2 in the 
5' end of the large subunit ribosomal gene (28S). 
PCR amplifications were performed in 25 /d 
reaction volumes and contained to a final concen¬ 
tration: lx Taq polymerase buffer, 2.5 mM MgClz, 
0.2 //M each primer, 0.8 mM dNTPs and 0.75U 
of Taq polymerase. The use of the hot start 
polymerase HotMaster Taq (Eppendorf) required 
an initial denaturation at 94°C for 2 min prior to 
the commencement of the remaining cycle param¬ 
eters; then followed 35 cycles of 94°C for 20 sec, 
55-58°C for 20 sec, 65°C for 45 sec and a final 
extension 65°C for 5 min. 

PCR products were gel purified using 'Perfect 
Prep' gel cleanup kit (Eppendorf) and forward 
and reverse sequencing reactions were carried 
out according to standard AB1 PRISM dye-deoxy 
terminator sequencing protocols using Big Dye 
Terminator versions 1.1 and 3.1. Chromatographs 
were checked and all sequences were aligned 
using Se-Al v2.0a!0 (Rambaut 1996). Estimates 
of sequence divergence including insertions 
(uncorrected p-distances) were calculated using 
the pairwise base distance function in PAUP* 
v4.0bl0 (Swofford 2002). We verified the origin 
of the amplified sequence data by conducting a 
BLAST search in GenBank thus determining 
the phylogenetic affinity with sequences from 
other actiniarian or anthozoan species. Sequences 
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for this same region were also obtained from 
GenBank from two individuals of the species 
Hetcmctis magnified (Accession no: AF050201 (H. 
magnified 1) and AF050211 (H. magnified 2)). 

SEA ANEMONE PRIMERS 

Based on the sequence results obtained from 
four of the six study species using the above 
described "universal' primers and one of the H. 
magnified sequences plus contaminating zooxan- 
thellae sequences from the remaining two species 
(Hetevaetis main and Maevodaetyla doveensis — 
see Results), we designed two new primers. 
These primers were designed to be specific to 
sea anemones and located in regions of identical 
sequence among the sea anemone species (for 
which we had data) but mismatched the zoo- 
xanthellae sequences at 45-50% of sites (see 
FIG. 1). These two new primers seaanemlSS : 5' 
TTA GTG AGG ACT CCT GAT TGG C 3' and 
seaanemlSS: 5' AGT CTC GCC TGA TCT GAG 
G 3' lie within 50bp downstream from RA2 and 
1 TS2.2 respectively. We tested the primers against 
the same six species used with the 'universal' 
primers. Amplification conditions, clean up and 
sequencing reactions with the new primers are 
identical to those described earlier. 

RESULTS 

DNA EXTRACTION 

In contrast to Pinto el al. (2000) no problems 
were experienced extracting DNA from ethanol 
preserved sea anemone tissues using the DNeasy 
tissue kit. Prior treatment of the samples to 
remove ethanol was not required; nor did the 
tissues need to be homogenised in liquid nitro¬ 
gen prior to the extraction process. Further¬ 
more, total tissue digestion was completed within 
1-3 hours at 55°C as recommended by the man¬ 
ufacture's protocol as opposed to the 72 hour 
period at 37»C used by Pinto et al. (2000). 
UNIVERSAL PRIMERS 

An 800bp (approximately) PCR fragment was 
successfully amplified from all six sea anemone 
species and all tissue types using the universal 
primers RA2 and ITS2.2. Readable sequence 
data of the fragment (including the 3' end of the 
18S gene, full length ITS1, 5.8S gene and 1TS2 
and the 5' end of 28S gene) was obtained from 
only three of the six species (?Aneinonia sp, M. 
doreensis and S. liaddoni). Partial/non-overlapping 


sequences were obtained from the remaining 
three species (E. quadricolor, H. main and O. 
ntuscosa). Not all tissue types generated readable 
sequence data. For example, sequences obtained 
from the pedal disc tissues of H. main and M. 
doreensis were unreadable with evidence of mul¬ 
tiple sequences present in the chromatograph 
(Table 2). This result was unexpected given that 
the amplified PCR product revealed a clear single 
band. However, readable sequence data was 
obtained from the tentacles of those same two 
species. BLAST searches of all readable sequences 
(either complete or partial) revealed strong match¬ 
es (90-97% identity) with other sea anemone 
and/or anthozoan species in GenBank for only 
four of the six study species (Table 2). The 
sequence data obtained from the tentacles from 
H. main and M. doreensis however, matched with 
almost 99% identity to other symbiotic algae 
sequences (e.g. Symbiodinium sp.) indicating prefer¬ 
ential amplification of the zooxanthellae DNA 
in each of these species. Interestingly, the sequence 
data obtained from both the pedal disc and 
tentacles of S. liaddoni were identical and BLAST 
searches of these and that obtained from the 
tentacles of E. quadricolor revealed closest simil¬ 
arity to other anthozoan species indicating that 
the host DNA had preferentially amplified and/ 
or that zooxanthellae are either not present or in 
high enough density to mask the host DNA in 
both these species. 

SEA ANEMONE PRIMERS 

Amplification success using our primers 
seaaucinlSS and seaaucm28S varied from that 
seen with the universal primers. Approximately 
750 bp were obtained from five of the six anemone 
species; no PCR product amplified from H. main 
regardless of tissue source (Fable 2). For the three 
species for which either tentacle and/or pedal 
disc tissues were available, amplification success 
varied from species to species. No PCR product 
was obtained from M. doreensis tentacle DNA; 
in contrast, product amplified from the tentacle 
DNA of E. quadricolor and both tissue types for 
S. liaddoni (Table 2). 

The lack of amplification success for H. main 
was surprising given that sequence data from 
the congeneric species, H. niagnifica, was used 
in the alignment from which the new primers 
were designed and that the regions of both the 
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Table 2. PCR and sequence results obtained from anemone tissues using both the 'universal' primers and sea 
anemone specific rDNA ITS primers. Presence (+) or absence (-) of product is indic ated . 




Preferential amplification and 

Preferential amplification and 



sequence obtained using 

sequence obtained using sea 



'Universal' primers 

anemone specific primers 

Species 

Tissue used in 

Anemone 

Zooxanthellae 

Anemone 

Zooxanthellae 

extractions 

DNA 

DNA 

DNA 

DNA 

lAnemonia sp. 

Column / 

Pedal disc 

+ 

- 

+ 

- 

Entacmaea 
\ quadricolor 

Tentacle 

+ 

- 

+ 

- 

Heteractls main 

Pedal disc 

+ 

+ 

- 

- 

Tentacle 

- 

+ 


— 

Macrodactyla 

Pedal disc 

+ 

+ 

+ 

- 

doreensis 

Tentacle 

- 

+ 

- 

- 

Oulactis mnscosa 

Column / 

Pedal disc 

+ 

- 

+ 

- 

Stichodactyla 

Pedal disc 

+ 

- 

+ 

- 

haddoni 

Tentacle 

+ 

- 

4 - 

- 


18S and 28S genes where these primers are located 
are identical among all the actiniarian genera 
(bar one site in O. wuscosa), for which sequence 
data was available. In order to see if we could 
amplify a product for H. walu but avoid zooxan- 
thellae DNA contamination, we tried the sea 
anemone primers in combination with the previ¬ 
ously successful universal primers; using seaanem 
18S paired with 1TS2.2 and scaanewISS paired 
with RA2. Successful amplification from H. wain 
DNA from both pedal disc and tentacles was 
only obtained using RA2/seaanew28S. 

Sequences, either partial or complete, obtained 
from lAnemonia sp., O. mnscosa and S. haddoni 
using the new sea anemone primers were iden¬ 
tical to those obtained using the universal primers, 
which had previously been confirmed as origin¬ 
ating from host anemone DNA rather than their 
algal symbionts. BLAST searches of complete 
sequences from M. doreetisis and H. wain obtained 
using anemone specific primers indicated greatest 
similarity to other anemones. Hence the anem¬ 
one specific primers had been successful in circum¬ 
venting the problems of zooxanthellae contam¬ 
ination. Curiously, E. qnadricolor did not return 
readable sequence data suggesting that further 
optimisation of the sequencing reaction for this 
species and these primers may be required. For 


H. wain , sequence obtained with RA2 revealed 
no mismatches in the 3' region of the 18S rRNA 
gene where seaanemlSS is located that would 
explain why this primer did not work on this 
species. Further experiments may be required 
to secure successful amplification with both 
anemone specific primers on this species. 

In summary, complete or overlapping sequen¬ 
ces of the 18S-28S fragment were obtained from 
only 3 species (lAnemonia sp. (724 bp), H. wain 
(670 bp) and S. haddoni (734 bp)). Although 
partial or non-overlapping sequences were ob¬ 
tained from E. qnadricolor (480bp from 18S), M. 
doreensis (523 bp: 300bp from 18S and 223bp 
from 28S) and O. wuscosa (556bp: 285bp from 
18s and 271bp from 28S), they were excluded 
from subsequent analysis due to incompleteness. 
SPECIES IDENTIFICATION 

Among the three species for which full 
sequences were obtained (including the two H. 
magnified sequences obtained from GenBank) 
estimates of sequence divergence ranged from 
0.14% within H. magnified up to 25.10% between 
H. magnified 1 and lAncwonia sp. (Table 3). The 
average level of sequence divergence among 
species was 23.84% indicating that this region 
may indeed prove to be useful for species 
identification in sea anemones. The exception 
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was the comparison between H. magnified and 
S. haddoni, where the divergence averaged only 
1.7% (Table 3). This result was somewhat unantici¬ 
pated given that it is significantly lower than the 
level of divergence found among the congeneric 
R magnified and R main sequences (ave 23.92%) 
and is therefore suggestive of possible taxonomic 
misidentifications. Considering thatR magnified 
or S. haddoni cannot be taxonomically verified 
for this study due to specimen/tissue unavail¬ 
ability, it highlights the importance of using 
genetics in conjunction with traditional taxonomic 
methods. 

The potential utility of this region for species 
identification is also evident from the example 
of lAnemonia sp., which was tentatively ident¬ 
ified in the field as resembling a juvenile 
Entacmaea quadrieolor collected from Bird Island 
(Table 1). While only partial sequences were 
obtained from the adult £. quadrieolor collected 
off Stradbroke Island, comparison of the 
sequences between the two specimens clearly 
showed they were significantly different (approx. 
18% sequence divergence over 480bp) and pos¬ 
sibly therefore two different species. Later histo¬ 
logical analysis revealed that the the Bird Island 
specimen was not £. quadrieolor as originally 
identified but may be TAnemonia sp., although 
the exact identity of this species still awaits final 
taxonomic confirmation. 

DISCUSSION 

The ribosomal DNA gene complex has proved 
highly successful for species identification across 
an incredibly broad range of taxonomic groups 
including plants (Chase el ai 2005), fungi 
(Kistaino et al. 1998; I wen el ai 2002), digencan 
parasites (Nolan & Cribb 2005) and mosquitos 


(Collins & Paskewitz 1996). It has even been 
used recently to identify commercial crustacean 
species from larvae collected in plankton sur¬ 
veys (Wang et ai 2006). In this study we inves¬ 
tigated for the first time, the utility of this region 
for identification of sea anemone species and 
the potential problems of using universal primers 
in species, which contain algal symbionts. 

While of a preliminary nature, our results 
showed high levels of sequence divergence 
among species using this region compared with 
divergence estimates an order of magnitude 
lower within a species indicating that it may 
indeed be ideal for assisting with sea anemone 
species identification. The questions at what 
taxonomic level and how useful this region may 
be for resolving phylogenetic relationships among 
sea anemone species was not the focus of this 
study but should certainly be investigated as 
more sequences become available. Acuna et ai 
(2007) used phylogenetic tools rather than 
estimates of sequence divergence to distinguish 
between different Aidactinia species and found 
extremely short branch lengths among individ¬ 
uals within a species compared to those between 
species. 

The usefulness of conserved 'universal' 
primers clearly depends on the species and 
tissue type available for analysis. However, as 
shown by the results obtained from R mailt and 
M. doreensis, extraction of 'uncontaminated' 
host DNA from samples taken only from pedal 
disc tissues clearly should never be assumed. In 
order to guarantee that host DNA is amplified 
alone, use of primers specific to sea anemones 
are recommended; if not on their own then at 
least in combination with another universal 
primer. The extent to which the primers designed 


Table 3. Estimates of sequence divergence among species for which the complete 18S-28S fragment was 
obtained (max 758bp). Sequences for H. magnified obtained from GenBank. * Specimen collected subsequent 
to Workshop. _ 



lAnemonia sp. 

II. magnified 1 

II magnified 2 

H. main 

S. haddoni 

lAnemonia sp. 

- 





H. magnified 1 

25.10% 

_ 




H. magnified 2 

24.93% 

0.14% 

- 



H. wain* 

21.60% 

23.84% 

23.99% 

~ 


S. haddoni 

24.42% 

1.64% 

1.77% 

22.98% 

- 
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for this study will work across all actiniarians 
remains to be seem Further preliminary PCR testing 
using seaane?ul8S and seacwem28S on another 
seven species from Moreton Bay, and repre¬ 
senting another seven actiniarian genera, proved 
highly successful with strong amplicons pro¬ 
duced in all seven species. Only subsequent 
sequencing will confirm whether or not the 
host DNA has been successfully targeted. 

Modem DNA extraction kits also seem highly 
useful for overcoming any difficulties associated 
with DNA extraction from ethanol preserved 
sea anemone tissues. Why we experienced so 
few problems compared with the earlier work 
of Pinto et at . (2000) is unclear. It may be that we 
were able to work with tissues from recently 
ethanol preserved specimens, rather than ones, 
that had been in ethanol for an extended time. 

Finally, a number of aspects of this study 
reinforce the value of being able to combine 
histological analysis with genetic testing to 
irrefutably verify a species' identity, especially 
given the embryonic stage of developing genetic 
markers for this group. In the case of a supposed 
juvenile E. quadricolor, the genetic data strongly 
indicated an incorrect field identification, and a 
subsequent histological analysis proved this to 
be so, identifying it instead as a probable Anemonia 
species. Furthermore, the curious result showing 
much greater sequence divergence between the 
two Heteractis species than that detected between 
Hetcractis magnified and Stichodactyla haddoni can¬ 
not, frustratingly, be resolved further. While 
again indicative of possible misidentifications, 
the H. magnified sequences available on GenBank 
are not associated with registered specimens 
and the S. haddoni cannot be analysed histo¬ 
logically as the whole animal was not collected 
from the field. 
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